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1.  HALL  COEFFICIENTS  IN  NARROWGAP  SEMICONDUCTORS 


We  have  completed  our  calculations  of  the  Hall  coefficients  in  HgCdTe  and  re¬ 
lated  materials.  The  Hall  fa^ctor  in  small  gap  semiconductors  is  calculated  by  solving 
Boltzmann  transport  equation  with  full  bandstructures  and  Fermi-Dirac  statistics. 
Previous  calculations  based  on  relaxation  time  approximations  indicated  a  large  de¬ 
viation  in  the  Hall  factor  from  unity.  However,  our  detailed  calculations  yield  a 
uniform  value  of  unity  for  all  temperatures  at  low  carrier  concentrations.  Surpris¬ 
ingly,  this  is  in  agreement  with  the  the  value  currently  used  in  experiments  which  is 
based  on  parabolic  approximation,  Maxwell- Boltzmann  distribution  and  relaxation 
time  approximation.  Further  details  of  this  work  are  given  as  Appendix  A. 

2.  DEFECT  MODELING  STUDIES  IN  HgCdTe  AND  CdTe 

An  invited  paper  entitled  “Defect  modeling  studies  in  HgCdTe  and  CdTe”  was 
presented  at  the  1994  MCT  Workshop  in  October.  This  paper  is  included  as  Ap¬ 
pendix  B. 

We  have  speculated  that  the  tellurium  antisite  is  both  the  residual  donor  and  a 
prinicpal  contributor  to  Shockley-Reed  recombination.  Thus  devising  means  to  re¬ 
ducing  the  density  of  this  defect  is  likely  to  improve  devices.  We  have  shown  that  the 
most  likely  diffusion  path  for  the  tellurium  antisite  is  through  a  mercury  vacancy  as¬ 
sisted  process.  Therefore,  to  decrease  the  tellurium  antisite  concentrations  one  must 
first  devise  annealing  strategies  to  accomplish  this  end  prior  to  the  low-temperature 
mercury-saturated  anneal  customairly  used  to  elimate  the  mercury  vacancy.  A  pre¬ 
scription  to  accomplish  these  goals  has  been  formulated,  and  experiments  to  test  its 
effectiveness  are  suggested. 

Another  major  accomplishment  reported  on  this  period  is  a  native  point  defect 
study  of  CdTe.  Once  again  we  agree  with  the  limited  experimental  information  avail¬ 
able,  and  the  details  we  present  suggest  means  to  imporve  the  material  properties. 

3.  TEMPERATURE  DEPENDENCE  OF  THE  BAND  GAPS  IN  HgCdTe  AND 
OTHER  SEMICONDUCTORS 

A  paper  on  the  temperature  dependence  of  the  semiconductor  band  gaps  was 
presented  at  the  1994  MCT  Workshop.  This  paper  is  included  as  Appendix  C. 
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In  the  previous  report  we  described  our  full-bandstructure  cedculation  of  Hall  factor, 
r,  as  a  function  of  impurity  concentrations  and  temperature  in  the  relaxation  time  approx¬ 
imation.  The  obtained  results  showed  a  considerable  increase  in  r,  with  a  maximum  value 
of  1.3,  is  expected  at  low  concentrations  and  high  temperatures.  These  observations  had 
a  strong  effect  on  experimentally  deduced  impurity  or  intrinsic  carrier  concentrations  and 
carrier  mobilities  zis  the  experimentalists  use  a  value  of  1  for  r.  In  this  period,  we  carried 
out  the  detailed  calculation  of  r  by  obtaining  full  solution  to  Boltzmann  transport  equa¬ 
tion  (BTE)  in  the  presence  of  both  electric  and  magnetic  fields.  In  all  our  calculations, 
accurately  obtained  hyperbolic  bandstructure  for  the  conduction  band  is  used. 

Several  methods  to  solve  the  BTE  with  full  bandstructures  and  Fermi-Dirac  (FD) 
statistics  are  described  in  the  literature.  However,  all  of  them  are  time  consuming  iterative 
schemes.  We  had  recently  extended  our  efficient  eigenvalue  approach  to  solve  BTE  with  FD 
statistics  and  external  electric  field^.  This  method  was  shown  to  be  fast  and  to  accurately 
obtain  the  electron  mobility  curves  for  various  impurity  concentrations  and  temperatures. 
We  further  extend  this  method  to  include  the  effect  of  crossed  magnetic  field. 

The  steady  state  BTE  in  the  presence  of  electric  (E)  and  magnetic  (B)  fields  is, 

V/(k) .  i{E  +  V  X  B)  =  2  (u.(k,k')/(k'){l  -  /(k))  -  Kik,  k)/(k)(l  -  /(k))]  (1) 

k' 

The  first  term  of  the  right  hand  side  is  the  gain  term  and  the  second  one  is  the  loss  term. 
In  equilibrium,  the  left  hand  side  of  Eq.(l)  is  identically  zero  and  because,  in  general  for  a 
system  of  interest  interacting  with  a  heat  bath  w(k,k’)  e"^®k'  =  w(k’,k)  e“^^it,  f  becomes 
the  equilibrium  FD  distribution  function  fo  given  by, 

/oo(£.)=(e‘^  +  l)’‘  (2) 

where  ep  is  the  Fermi  energy.  The  solution  f(k)  can  always  be  written  as  a  sum  of  foo(k) 
and  some  deviation,  A(k).  That  is, 

/(k)  =  /oo(k)-hA(k)  (3) 

Note  that  A(k)  is  zero  for  a  sample  with  ohmic  contacts  because  the  number  of  elec¬ 
trons  should  remain  constant.  Combining  Eq.(l)  through  Eq.(3),  we  get, 

V/(k)  •  ^(E  +  V  X  B)  =  I^[H'(k,  k)A(k)  -  IVCk,  k)A(k) 

+  (ti.(k,k)-u.(k,k))A(k)A(k)]  (4) 
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where  the  renormalized  and  the  usual  transition  probability  per  unit  time  w  are  related 

by, 


=  (5) 

Note  that  for  elastic  scattering  W  and  w  axe  equal.  However,  for  inelastic  cases  the 
relative  size  of  to  w  depends  on  whether  energies  at  k  and  k  are  larger  or  smaller 
than  cf  •  If  both  initial  and  final  energies  are  larger  (or  smaller)  than  cf  ,  then  only 
small  corrections  to  w  occur.  However,  if  the  initial  state  is  above  cf  and  the  final  state 
is  below  Cf  ,  then  for  that  scattering  event  W  is  suppressed.  This  tends,  for  example,  to 
decrease  the  contribution  of  inelastic  scattering  events  involving  phonon  emission.  Also 
notice  that  in  the  event  the  perturbation  is  small,  then  the  A*  term  can  be  neglected  and 
Eq.(4)  resembles  the  traditional  MB  steady  state  BTE  where  W  plays  the  role  of  w. 

We  further  expand  A(k)  in  power  series  of  an  indexing  parameter  A  and  /3,  which  we 
will  eventually  set  to  unity.  Thus, 


A(jfc)=  f;  /„m(k)A"^'"  (6) 

n,m=l 

Substituting  Eq.(6)  in  Eq.(4)  and  noting  that  f(k)  has  one  additional  term,  we  get. 


E  V/„,„(k).^(E  +  vxB)A"/3'"  = 

n,mssO 

f;  [H'Ck,  -  H'(k',  k)/„„(k)] 

+  [u;(k',k)-u;(k,k')]/„m(k)/„-,n-(k')A"+"'/3-+'"'  (7) 

Noting  that  the  indices  A  and  /?  are  proportional  to  E  and  B  respectively,  we  can  collect 
the  coefficients  of  the  same  power  of  A  and  /3  on  both  sides  to  get  a  series  of  equations. 
However,  for  our  purpose  of  calculating  r  where  only  small  fields  are  applied,  we  consider 
only  X°/3^,  and  A^/?^  terms.  We  get, 

V/oo(k)  •  =  5:  [W(k,  k')/io(k')  -  W(k',  k)/io(k)]  (8a) 

k' 

V/„(k)  4(v  X  B)  =  5;  [(^(k,  k')/o,(k')  -  ty(k',  k)/.,(k)]  (86) 

k' 
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and 


V/io(k)  ■  5(v  X  B)  =  2  [H'(k,  k')/i.{k')  -  H'Ck',  k)/„(k)]  (8c) 

k' 

In  the  isotropic  case  that  we  study  V/oo(k)  and  vare  proportional  to  kand  hence  from 
Eq.(8b)  foi  is  identically  zero.  Note  that  this  information  is  used  in  obtaining  Eq.(8c). 
Knowing  foo  (from  Eq.(2)),  Eq,(8a)  is  solved  for  fio  and  the  solution  is  used  Eq.(8c)  for 
fii.  This  procedure  can  be  continued  to  the  required  precision.  We  emphasize  that  for  low 
field  electric  fields,  only  the  lowest  order  equations  given  here  need  to  be  solved  to  obtain 
accurate  answers.  When  hot  electron  effects  axe  addressed,  iterative  solutions  should  be 
used. 

As  shown  in  previous  publication  (Krishnamurthy  and  Sher,  J.  Appl.  Phys.,  75,  7904 
(1994)),  Eqs.(8a)  and  (8c)  are  efficiently  solved  in  a  matrix  inversion  method  by  expanding 
the  distribution  functions  in  terms  a  finite  set  of  basis  functions.  Once  the  functions  fio 
and  fii  are  obtained,  the  drift  velocity,  v<{,  drift  mobility,  /i<jj  Hall  coefficient,  R  and  Hail 
factor,  r  are  easily  calculated.  We  assume  that  E  is  in  z  direction  and  B  is  in  y  direction. 


ne 


<7  =  nefiE 
Je  =  nevf 
js  =  nevf 


.,E  _ 

Ek'/oo(k') 

,,B  _  Ek^k/u(^) 

Ek'/oo(k') 


All  of  the  above  equations  can  be  combined  to  yield  a  final  expression. 
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r=-ii- 

HEvfB 

We  calculated  the  Hall  factor,  r  for  GaAs  and  Hgo.TsCdo.jjTe  alloy  as  a  function  of  tem¬ 
perature  and  impurity  concentrations.  Interestingly,  the  calculated  values  at  low  carrier 
concentrations  are  always  around  1  for  any  temperature.  While  the  final  result  comfort¬ 
ably  justifies  the  value  currently  used  by  experimentalists,  it  should  be  noted  that  physics 
is  quite  different.  The  value  of  1  for  r  is  obtained  with  parabolic  bandstructures,  Maxwell- 
Boltzmann  distribution  and  relaxation  time  approximation.  We  had  shown  that  all  these 
approximations  are  incorrect  for  small  gap  semiconductors  like  HgCdTe.  While  the  re¬ 
moval  of  approximations  correct  the  mobilities,  Fermi  energy  and  temperature  dependent 
band  gap  substantially,  the  Hall  factor  remain  unaffected. 

We  are  in  the  process  of  calculating  Hall  factor  for  high  dopant  concentrations. 
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DEFECT  MODELESTG  STUDIES  IN  HgCdTe  and  CdTe 


M.  A.  HERDING,  A.  SHER,  and  M.  VAN  SCHILFGAARDE 

SRI  International,  Menlo  Park,  California  94025 

ABSTRACT 

We  have  used  a  quasichemical  formalism  to  calculate  the  native  point  defect  densities 
in  X  =  0.22  Hgi_,Cd*Te  and  CdTe.  The  linearized  muffin-tin  orbital  method,  based 
on  the  local  density  approximation  and  including  gradient  corrections,  has  been  used  to 
calculate  the  electronic  contribution  to  the  defect  reaction  free  energies,  and  a  valence 
force  field  model  has  been  used  to  calculate  the  changes  to  the  vibration  free  energy  when 
a  defect  is  created.  We  find  the  double  acceptor  mercury  vacancy  is  the  dominant  def^t, 
in  agreement  with  previous  interpretations  of  experiments.  The  tellurium  antisite,  which 
is  a  donor,  is  also  found  to  be  an  important  defect  in  this  material.  The  mercury  vacancy 
tellurium  antisite  pair  is  predicted  to  be  well  bound,  ^d  is  expected  to  be  important  for 
tellurium  antisite  diffusion.  We  consider  the  possibilities  that  the  tellurium  antisite  is  the 
residual  donor  and  a  Shockley-Read  recombination  center  in  HgCdTe,  and  sugg^tions 
for  further  experimental  work  are  meide.  We  predict  that  the  cadmium  interstitial,  a 
double  acceptor,  is  the  dominant  defect  for  low  cadmum  pressures,  while  the  cadmium 
interstitial,  a  double  donor,  dominates  at  high  cadmium  pressures. 

Key  words:  HgCdTe,  CdTe,  defects,  defect  complexes 

INTRODUCTION 

The  pseudobinary  semiconductor  alloy  Hgi_*Cdj.Te  with  x=0.22  is  currently  the  ma¬ 
terial  of  choice  for  high-performance  detectors  in  the  long-wavelength  infrared  (LWIR) 
(8-14  ^m).  Unlike  most  other  II-VI  systems,  both  extrinsic  p-  and  n-type  doping  can 
be  achieved  in  HgCdTe,  although  in  as-grown  material  the  electrical  characteristics  are 
often  determined  by  native  point  defect  concentrations.  Understanding  the  properties  of 
point  defects  and  manipulation  of  their  densities  during  growth  aud  processing  is  essential 
to  high-yield  manufacturing  of  focal  plane  arrays  (FPAs).  As  in  other  semiconductors, 
it  is  difficult  to  establish  the  presence  and  identity  of  all  the  important  neutral  and  com¬ 
pensating  point  defects  during  growth  and  processing,  much  less  to  determine  their  con¬ 
centrations.  CdTe  is  important  both  as  a  substrate  and  passivating  material  for  epitaxial 
layers  of  LWIR  HgCdTe.  Native  point  defects  are  of  interest  in  CdTe  in  that  they  relate  to 
its  stoichiometry  (which  in  turn  hzis  been  shown  to  impact  the  minority  carrier  lifetimes^ 
in  HgCdTe  devices),  the  formation  and  annihilation  of  tellurium  precipitates and  the 
self—  and  inter— diffusion  coefficients  that  impact  materials  stability  during  growth,  during 
subsequent  processing,  and  over  the  device  lifetime. 

Our  goal  in  this  paper  is  to  theoretically  identify  the  important  native  defects  in 
HgCdTe  and  CdTe,  to  calculate  their  densities  as  a  function  of  growth  and  processing 
conditions,  to  validate  our  predictions  by  comparison  with  experimentally  deduced  prop¬ 
erties  of  the  native  defects,  and  to  suggest  new  experiments  to  begin  to  unravel  the 
remaining  mysteries  in  these  materials.  We  have  included  in  our  analysis  of  HgCdTe  both 
neutrcd  and  ionized  states  of  eight  native  point  defects  and  one  defect  pair.  Our  focus  is  on 
X  =  0.2  (for  comparison  to  annealing  data®),  x  =  0.22  (for  LWIR  applications),  and  x  =  1 
(for  substrate  and  passivating  layers)  Hgi_*Cd,Te.  Our  analysis  of  CdTe  includes  eight 
neutral  native  point  defects.  As  we  will  show,  we  have  attempted  to  incorporate  all  of  the 
important  contributions  to  the  defect  formation  free  energies  and  adopt  a  first— principles 
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Approach  for  inost  of  the  quantities  we  calculate. 

METHODOLOGY 

To  calculate  the  concentration  of  the  native  point  defects  in  HgCdTe  we  employ  a 
quasichemical  formalism'*  in  which  the  formation  of  independent  crystalline  defects  is 
expressed  as  chemical  reactions.  For  example,  for  the  formation  of  the  neutral  mercury 
vacancy  in  HgTe  with  the  mercury  vapor  as  the  reference  phase,  we  have  the  reaction 

HgTe-.Virje  +  Hg„,,^.  (1) 

Our  notation  is  that  of  Kroger^  in  which  the  primary  symbol  is  the  defect  species  and 
the  subscript  indicates  the  site  the  defect  occupies;  V  indicates  the  vacancy,  and  /  the 
interstitial,  and  no  subscript  indicates  the  species  is  occupying  its  normal  lattice  position; 
X  indicates  the  neutral  defect  species.  Although  we  have  chosen  the  mercury  in  the  vapor 
phase  as  our  reference  and  thereby  will  choose  a  mercury  pressure  to  determine  the  chem¬ 
ical  potential  of  mercury  in  the  system,  one  could  have  chosen  some  other  reference  state 
for  the  mercury  or  tellurium  (for  example,  by  specifying  a  binary  solution  of  Hgi-^Te,  at 
some  temperature  T  in  coexistence  with  the  HgTe  solid,  as  is  essentially  done  during  liq¬ 
uid  phase  epitaxy  (LPE)  growth).  For  our  equilibrium  calculations,  the  lirmts  of  mercury 
pressures  within  the  existence  region  of  Hgi_rCd*Te  are  taken  from  experunent.® 

In  a  similar  manner,  reactions  can  be  written  for  the  other  point  defects  of  the  system. 
In  our  analysis  we  have  included  eight  native  point  defects  (plus  their  ionized  species 
and  the  electron  and  hole):  the  mercury  and  tellurium  vacancies  (Vag  and  Vt«),  the 
mercury  and  tellurium  antisites  (Hgx*  aJid  Teag,  and  two  types  of  mercury  and  tellurium 
tetrahedral  interstitials — one  surrounded  by  four  cation  near— neighbors  (Hg/„j  and  Te/„,) 
and  the  other  surrounded  by  four  anion  near-neighbors  (Hg/^.  and  Te/^.)-  We  have  also 
included  the  bound  mercury  vacancy  tellurium  2uitisite  pair  ( Vng  —  T^g)pair  y  in  which  the 
vacancy  eind  antisite  occupy  near-neighbor  cations  sites,  via  the  reaction 

Vng  +  TeHg  (Vkg  -f  TeHg)p^j^ .  (2) 

For  low  densities  of  noninteracting  defects,  the  law  of  mass  action  can  be  used  to 
determine  the  neutral  defect  concentrations.  For  the  mercury  vacancy  this  corresponds 
to 


where  Fyx  is  the  reaction  free  energy  for  the  neutral  mercury  vacancy,  ks  is  Boltzmann’s 
constant,"r  is  the  temperature  in  kelvin,  and  0  is  the  number  of  unit  cells  per  unit  volume. 
Thus,  once  the  reaction  free  energy  is  known,  it  is  straightforward  to  calculate  the  defect 
concentrations.  Of  course,  all  the  work  is  involved  in  the  calculation  of  the  reaction  free 
energies. 

We  have  attempted  to  calculate  all  of  the  important  contributions  to  the  defect  reac¬ 
tion  free  energies.  An  electronic  contribution  to  the  free  energy  results  from  the  change 
in  the  total  electron  energy  of  the  sohd  when  a  neutral  defect  is  created;  included  in 
this  energy  is  the  electronic  energy  of  free  mercury  atoms  generated  or  consumed  in  the 
defect  reaction.  To  calculate  the  electronic  contribution  we  employ  the  self-consistent 
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first-principles  full-potential  linearized  muffin-tin  orbital  method®  and  the  local-density 
approximation  (LDA),  including  gradient  corrections  of  the  Langreth-Mehl-Hu  type’^ 
(which  greatly  improves  the  overbinding  foimd  in  the  LDA*-®).  The  vibrational  modes 
of  the  system  are  also  changed  when  a  defect  is  created;  we  calculate  this  change  to  the 
formation  free  energy  (both  the  enthalpy  and  entropy)  using  a  Green’s  function  formal¬ 
ism  within  a  valence  force  model  plus  point-charge  ionic  model.  An  entropy  contribution 
to  the  formation  free  energy  also  arises  from  the  partial  occupation  of  degenerate  lev¬ 
els  associated  with  the  defect  and  from  the  introduction  of  a  preferred  direction  via  a 
symmetry-lowering  Jahn-Teller  distortion.^®  The  combination  of  the  electronic,  vibra¬ 
tional,  degeneraicy,  and  translational  (for  the  calculation  of  the  chemical  potential  of  the 
mercury  atom  in  the  vapor  phase)  free  energies  encompasses  the  primary  contributions 
to  the  total  defect  formation  free  energies  when  referenced  to  a  mercury  vapor.  Details  of 
the  results  for  these  energies  are  given  in  Ref.  9.  We  have  completed  only  preliminary  cal¬ 
culations  for  the  binding  energy  of  the  mercury  vacancy  tellurium  antisite  pair  indicated 
in  Eq.  3,  using  a  32-atom  supercell  and  with  overall  lattice  constant  relaxation  only.  We 
expect  that  this  defect  may  show  further  relaxation,  which  could  lower  the  defect  pair 
binding  energy  and  therefore  increase  their  density. 

Although  we  are  interested  in  studying  the  properties  of  point  defects  in  HgCdTe, 
the  neutral  defect  reaction  free  energies  are  calculated  for  HgTe;  because  we  predict  that 
the  defects  with  highest  concentrations  (the  mercury  vacancy  and  the  tellurium  antisite) 
reside  on  the  cation  sublattice  and  therefore  have  four  tellurium  atoms  for  first  nearest 
neighbors,  we  expect  that  the  presence  of  ceidmium  will  introduce  a  minor  modification 
to  the  formation  free  energies,  and  of  the  order  of  other  approximations  made  in  the 
calculation  (for  exeimple,  using  superceUs).  The  presence  of  cadmium  may  have  a  larger 
impact  on  interstitial  formation  free  energies,  which  can  occupy  sites  with  cations  as  near¬ 
neighbors.  Although  we  find  the  density  of  interstitials  to  be  relatively  low,®  the  cation 
interstitial  in  particular  is  very  important  in  annealing  of  HgCdTe.**  For  the  present  we 
have  included  the  effects  of  the  cadmium  only  by  adjusting  the  number  of  sites  available 
for  mercury  vac£incies  and  in  calculating  of  the  band  structure  used  in  predicting  the 
ionized  point  defect  concentrations.  For  further  discussion  and  details  of  the  Ccilculation 
of  the  reaction  free  energies,  see  Ref.  9. 

Because  the  native  point  defects  wiU  in  general  have  localized  levels  in  the  band 
gap,  we  need  to  calculate  the  concentration  of  the  ionized  defects  in  addition  to  the 
neutral  defects  discussed  above.  Once  the  density  of  the  neutral  defects  is  determined, 
the  concentration  of  ionized  defects  can  be  calculated  from 
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for  a  donor  defect  Y.  Once  again  we  have  adopted  the  notation  of  Kroger:  the  bullet 
superscript  indicates  a  positive  charge  and  a  prime  a  negative  charge.  Ea  and  Ed  are 
the  one-electron  acceptor  and  donor  levels;  fiF  is  the  Fermi  energy,  which  is  determined 
by  demanding  charge  neutrality;  2ind  gx  is  the  degeneracy  of  the  state,  including  both 
spin  and  Jahn-Teller*®  splitting.  Generalizations  of  these  expressions  are  used  for  multiply 
ionized  defects.  The  calculation  of  the  electron  and  hole  populations  demands  a  knowledge 
of  the  temperature-dependent  band  structure  and  use  of  Fermi-Dirac  statistics.  The  only 
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significiint  empirical  data  we  employ  in  this  calculation  are  those  needed  to  obtain  the 
temperature-dependent  band  structure,  which  is  extrapolated  to  high  temperatures  at 
which  equilibration  and  annealing  take  place  and  for  which  we  are  going  to  predict  defect 
concentrations.  Further  details  are  in  Ref.  9. 

DEFECTS  IN  HgCdTe 

Equilibrium  native  point  defect  concentrations:  annealing  studies 

Gibbs’  phase  rule  tells  us  that  for  a  system  of  three  components  (in  our  case  mercury, 
cadmium,  and  tellurium)  and  two  phases  (zincblende  solid  and  vapor)  there  are  three 
degrees  of  freedom.  In  evaluating  the  equilibrium  defect  concentrations  in  Hgo  sCdo  jTe 
we  have  chosen  the  temperature,  the  mercury  pressure  Pug,  and  the  alloy  composition  x  as 
these  specified  variables;  the  tellurium  and  cadmium  pressures,  the  crystal  stoichiometry, 
and  the  density  of  the  various  native  point  defects  are  determined  by  these  parameters. 

In  Fig.  1  we  show  our  results  for  the  77  K  hole  concentrations  in  x  =  0.2  Hgi_,Cd,Te  as 
a  fimction  of  mercury  partial  pressure  for  various  high— temperature  annealing  conditions, 
and  compare  them  to  the  results  of  Vydyanath.^-^®  We  have  assumed  that  the  high- 
temperature  defect  structure  is  frozen  in  during  a  quench  to  77  K.  All  eight  point  defects 
discussed  above  have  been  included  in  our  analysis.  The  agreement  with  experiment  is 
quite  good,  given  that  the  only  empirical  data  used  were  of  the  temperature  dependence 
of  the  band  structure  used  to  calculate  the  intrinsic  reaction  constant.  To  demonstrate 
the  sensitivity  of  our  results  to  the  free  energies  we  are  calculating  from  first  principles,  we 
have  cilso  shown  in  Fig.  1  the  results  of  a  calculation  using  a  mercury  vacancy  formation 
energy  that  has  been  increased  roughly  10%,  plus  a  constant  multiphcative  constant  of 
the  low— temperature  hole  concentrations  at  all  pressures  and  temperatures;  as  one  can 
see,  such  minor  modifications  to  our  calculated  parameters  yield  low— temperature  hole 
concentrations  that  are  in  very  good  agreement  with  experiment.  Certainly,  refinements 
to  our  theory  (for  example,  including  a  more  accurate  high-temperature  band  structure, 
more  precise  incorporation  of  alloy  effects  on  formation  free  energies,  anharmonic  effects 
in  the  vibrational  free  energies,  and  going  beyond  the  local  density  approximation)  could 
account  for  discrepancies  of  this  magnitude.  In  addition,  the  experimental  data  may  be 
impacted  by  reequihbration  during  the  quench  from  high  temperature. 

At  all  temperatures,  we  predicted  that  equilibrated  material  will  be  p— type  and  that 
the  dominant  defect  is  the  doubly  ionized  mercury  vcicancy,  in  agreement  with  mobility 
data.®  Our  predictions  differ  from  the  analysis  in  Ref.  3  in  several  ways,  though.  First, 
although  the  data  indicate  that  the  material  is  intrinsic  at  all  annealing  temperatures, 
we  predict  that  the  material  will  be  extrinsic  for  the  higher  anneahng  temperatures  with 
[h*]  ~  while  at  the  lower  annealing  temperatures  the  material  will  be  intrinsic  with 

[h*]  ~  (for  details  see  Ref.  9).  This  may  be  due  to  either  a  vacancy  formation  energy 
that  is  too  small  in  our  calculation  (as  demonstrated  by  the  dashed  fine  in  Fig.  1)  or 
errors  in  the  high-temperature  band  structure,  which  we  have  extrapolated  from  the  low- 
temperature  formulas.*'*’^®  Second,  we  find  at  higher  temperatures  and  lower  pressures, 
that  the  hole  concentrations  increase  with  increasing  pressure,  contrary  to  what  would  be 
expected  due  to  mercury  vacancy  acceptors  alone.  In  our  theoretical  predictions  this  is 
due  to  compensation  from  tellurium  antisites,  which  we  predict  are  donors,  and  which  are 
the  second  most  dominant  native  point  defect  in  Hgi_*CdxTe.  Because  the  equiUbrium 
cintisite  density  varies  roughly  as  compared  to  P^^  to  P^^j^  for  the  mercury  vacancy, 
it  is  most  important  at  low  mercury  pressures.  At  the  highest  temperatures  there  does 
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appear  to  be  a  corresponding  role-ofF  in  the  experimental  data  at  the  lower  pressures, 
which  may  be  indicative  of  compensation  by  tellurium  antisites,  but  it  could  also  be  an 
experimental  artifact  caused  by  quenching  inefficiencies  for  the  higher  temperatures. 

Low-temperature  mercury-saturated  anneals  are  of  technological  importance  to  re¬ 
ducing  the  as-grown  p-type  carrier  concentrations  or  to  convert  the  material  to  n-type  in 
nominally  imdoped  material.  In  Fig.  2  we  have  plotted  the  concentrations  of  mercury  va¬ 
cancies  and  tellurium  antisites  as  a  function  of  inverse  temperature  for  constant  pressures 
within  the  existence  region;  neither  of  these  concentrations  included  the  defects  that  are 
bound  into  mercury  vacancy  tellurium  antisite  pairs.  The  lower  boundaries  correspond 
to  the  defect  densities  for  mercury-saturated  anneals,  the  upper  boundary  to  tellurium- 
saturated  anneals.  As  one  can  see,  at  a  given  pressure  and  temperature,  the  mercury 
vacancy  concentrations  are  always  in  excess  of  the  tellurium  antisite  population,  and  in 
equilibrium  the  material  should  always  be  p-type.  As  discussed  in  the  following  section, 
we  expect  that  full  equilibration  of  the  tellurium  antisite  density  may  not  take  place  at 
lower  annealing  temperatures,  so  that  the  equilibrium  concentrations  predicted  for  this 
defect  may  not  always  be  realized. 

Annealing  schedules  are  also  of  technological  importance  in  forming  p-n  junctions  in 
as— grown  p— type  material  via  mercury  in  diffusion.^®'*^  The  results  of  these  experiments 
depend  on  the  mechanism  by  which  mercury  diffuses  and  are  related  to  the  identity  of 
the  residual  donor,  both  of  which  will  be  discussed  further  below. 

The  tellurium  antisite  in  Hgo.78Cdo.22Te 

Having  predicted  that  the  tellurium  antisite  will  be  present  in  substantial  densities, 
we  now  address  the  evidence  for  their  existence  in  LWIR  HgCdTe.  We  have  predicted 
that  the  tellurium  antisite  is  a  donor,  although  we  have  not  yet  extended  our  theory 
to  predict  whether  it  is  a  single  or  double  donor,  nor  the  precise  location  of  the  defect 
levels  in  the  gap.  This  is  perhaps  the  most  difficult  part  of  our  calculation  because  of 
the  LDA  band  gap  errors  intrinsic  to  our  method  and  the  dispersion  arising  from  the  use 
of  supercells  that  give  defect— level  widths  larger  than  the  Hgo.78Cdo.22Te  band  gap.  To 
establish  the  presence  of  tellurium  antisites  without  a  prediction  of  the  donor  level(s)  in 
the  gap,  we  must  consider  the  possibility  that  there  are  resonant  level(s),  shallow  level(s), 
deep  level(s),  or  some  combination  of  these,  associated  with  the  tellurium  antisite,  and 
look  for  evidence  supporting  the  presence  of  the  antisite. 

One  of  the  outstanding  questions  in  the  properties  of  HgCdTe  is  the  identification 
and  elimination  of  the  residual  donor  that  controls  the  carrier  concentration  in  material 
annealed  at  low  temperatures  under  mercury-saturated  conditions.  Although  purification 
of  starting  materials  has  led  to  a  lowering  of  the  residual  donor  density,  a  lower  limit  of 
~  10^'*  cm~®  seems  to  have  been  reached.  To  examine  the  possibihty  that  the  tellurium 
antisite  is  the  residual  donor,  we  have  calculated  the  defect  densities  at  500*C,  roughly  the 
LPE  growth  temperature  from  both  the  mercury-  and  tellurium-rich  melts,  at  185®C,  the 
typici  molecular  beam  epitaxy  (MBE)  growth  temperature,  and  at  220*  C,  roughly  where 
low— temperature  mercury— rich  anneals  are  done;  results  are  shown  in  Fig.  3(a-c).  It  is 
clear  that  the  material  is  expected  to  be  p-type  for  both  MBE-  and  LPE-  grown  material, 
and  that  although  postgrowth,  mercury-saturated  anneals  (the  right  side  of  Fig.  3(c)) 
will  lower  the  vacancy  density  below  the  10^^  cm“^  level,  the  material  is  still  predicted  to 
be  mercury— vacancy— doped  p— type.  We  thus  conclude  that  the  residual  donor  is  not  due 
to  an  equihbrium  concentration  of  native  point  defects. 

The  possibility  remains  that  the  residual  donor  is  associated  with  a  nonequilibrium 
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defect  concentration  of  tellurium  antisites.  We  expect  that  the  diffusion  coefficient  of 
the  tellurium  antisite  is  relatively  small  because  the  diffusion  of  an  antisite  will  neces- 
garily  involve  at  least  one  additional  point  defect,  such  as  the  mercury  vacancy  or  the 
tellurium  interstitial  (diffusion  via  a  consorted  exchange  of  a  tellurium  antisite  with  a 
mercury  atom  on  an  adjacent  cation  lattice  site  is  unlikely).  As  such,  the  tellurium  an¬ 
tisite  may  not  reach  equilibrium  densities  for  the  times  and  temperatures  corresponding 
to  the  low-temperature  mercury-saturated  anneals  typically  employed  to  equihbrate  the 
mercury  vacancy  density  (Fig.  3(c)).  If  tellurium  antisite  densities  are  in  fact  equilibrated 
at  a  temperature  at  which  the  antisite  diffusion  effectively  stops  during  cooldown  from 
the  growth  temperature,  then  the  antisites  may  be  frozen  in  at  higher,  nonequilibrium 
concentrations  corresponding  to  the  residual  donor  density. 

Addressing  the  question  as  to  why  the  antisite  density  might  be  frozen  in  at  roughly 
the  same  concentration  for  LPE  material  grown  from  both  the  mercury  and  tellurium 
melt,  we  return  to  the  means  by  which  tellurium  antisite  diffuses  in  the  material,  and 
assume  it  diffuses  by  a  mercury  vacancy  mechanism.  This  assumption  is  motivated  by 
several  factors.  First,  this  mechanism  involves  only  one  point  defect  in  addition  to  the 
tellurium  antisite,  and  involves  a  simple  migration  mechanism  between  the  two  defects. 
Second,  it  involves  the  mercury  vacancy,  whose  density  is  fairly  high,  and  therefore  the 
probabihty  of  tellurium  antisite  mercury  vacancy  pairs  is  expected  to  fairly  high.  In 
addition,  the  mercury  vacancy  is  an  acceptor  and  the  tellurium  antisite  is  a  donor  so 
they  are  expected  to  form  a  bound  pair  based  on  Coulombic  attraction,  and  the  mercury 
vacancy  is  too  small  for  the  lattice,  while  the  tellurium  antisite  is  too  large,  so  that  there  is 
a  mechanical  attraction  between  them  as  well.  Finally,  the  migration  mechanism  involving 
an  interstitial — for  example,  via  a  kickout  mechanism  whereby  a  mercury  interstitial  kicks 
out  the  tellurium  antisite  to  form  a  tellurium  interstitial — ^is  unlikely  since  it  involves 
defects  that  are  all  too  large  for  the  lattice  and  are  donors  and  therefore  are  unlikely 
to  form  pairs.  The  kickout  mechanism  may  be  important  when  mercury  interstitials 
are  injected  into  the  material,  and  is  discussed  later  in  this  paper.  The  likelihood  of 
the  tellurium  antisite  diffusion  proceeding  by  the  mercury  vacancy  mechanism  is  further 
supported  by  our  prehminary  prediction  of  the  mercury  vacancy  tellurium  antisite  binding 
energy  of  1.1  eV;  the  corresponding  densities  are  shown  in  Fig.  3.  We  show  our  predictions 
for  the  pair  density  as  a  range  in  which  the  lower  hmit  corresponds  to  the  defect  pair  being 
neutral,  as  our  preliminary  predictions  indicate,  and  the  upper  hmit  corresponding  to  the 
pair  having  a  donor  state  at  the  valence  band  edge  and  an  acceptor  level  at  the  conduction 
band  edge.  Unhke  the  native  point  defects,  we  expect  that  these  pair  defect  concentrations 
may  change  as  we  refine  the  free  energy  calculations.  If  tellurium  antisite  diffusion  is  via 
the  mercury  vacancy,  the  diffusion  coefficient  will  be  proportional  to  the  density  of  the 
mercury  vacancy  tellurium  antisite  pairs,  that  is 

^T«h,  ~  [(  Vhg  +  TeHg)patr]  •  (6) 

From  Fig.  3(a)  for  LPE  growth,  one  can  see  that  the  density  of  defect  pairs  for  low 
mercury  pressures,  corresponding  to  material  grown  from  the  tellurium-rich  melt,  is  four 
orders  of  magnitude  higher  than  for  that  grown  from  the  mercury-rich  melt.  This  implies 
that  although  material  grown  from  the  tellurium  melt  will  contain  higher  densities  of 
tellurium  antisites,  in  the  cooldown  from  the  growth  temperature  the  tellurium  antisites 
will  reequihbrate  much  more  rapidly  due  to  the  large  pair  density  th£m  will  material  grown 
from  the  mercury-rich  melt.  In  contrast,  for  material  grown  from  the  mercury  melt,  the 
tellurium  antisites  present  at  the  growth  temperature,  although  lower,  may  be  frozen  in 
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because  the  pair  density  is  also  lower. 

We  now  turn  to  growth  by  MBE,  which  takes  place  at  185*C  under  mercury-poor 
conditions,  corresponding  to  the  left  side  of  Fig.  3(b).  At  the  phase  stability  boundary, 
we  predict  that  the  material  will  be  mercury-vacancy  doped  p-type,  with  a  carrier  con¬ 
centration  of  ~  10^®  cm-®.  Experimentally  as-grown  material  is  found  to  be  either  n-  or 
p-type  with  carrier  concentrations  in  the  10*®  cm"®  range.*®  Because  MBE  growth  is  a 
nonequilibrium  process,  it  is  possible  that  growth  may  take  place  beyond  the  existence 
region;  for  example,  if  the  equivalent  pressure  were  of  the  order  of  10“^  atm  at  185®C,  the 
materials  would  be  tellurium  antisite  doped.  While  this  provides  a  possible  explanation 
of  how  the  material  could  be  n-  or  p-type  as  grown  based  on  equilibrium  concentrations 
of  defects,  extending  this  argument,  one  would  expect  under  some  growth  conditions  to 
be  able  to  obtain  highly  compensated  material  with  very  low  carrier  concentrations;  this 
is  never  seen,  to  our  knowledge.  A  more  likely  possibility  is  that  nonequilibrium  densities 
of  point  defects  are  frozen  into  the  MBE  material  resulting  from  details  of  the  surface 
kinetics. 

There  is  substantial  experimental  evidence  of  deep  levels  in  LWIR  HgCdTe,  which 
have  not  yet  been  definitively  associated  with  any  particular  defect.  The  nunority  carrier 
lifetimes  in  vacancy— doped  material  are  limited  by  Shockley— Read  recombination;  deep- 
level  transient  spectroscopy  (DLTS)  measurements*®-®®-®*  have  identified  two  donor-like 
defect  levels  at  0.4E,„p  and  OJbEgap-  The  densities  of  these  levels  roughly  track  the  mer¬ 
cury  vacancy  concentration,  although  substantial  scatter  in  the  correlation  is  observed. 
Neither  of  these  levels  is  believed  to  be  associated  with  the  mercury  vacancy  itself,  nor  do 
they  appear  to  be  associated  with  the  same  defect  because  they  do  not  track  one  another. 
In  equilibrium,  the  tellurium  antisite  density  will  track  with  the  mercury  vacancy  con^n- 
tration,  although  it  will  show  a  sharper  dependence  on  the  mercury  pressure.  In  addition, 
as  discussed  above,  we  do  not  expect  that  equilibrium  concentrations  of  tellurium  antisites 
will  be  present  except  at  very  high  temperatures  corresponding  to  liquid-phase  growth, 
and  thus  the  ratio  of  tellurium  antisites  to  mercury  vacancies  expected  in  equilibrium 
may  not  be  experimentally  realized,  leading  to  substantial  deviations  from  the  equilib¬ 
rium  ratio,  and  scattering  in  their  concentrations.  Thus,  it  is  plausible  that  one  of  the 
Shockley— Read  recombination  centers  is  associated  with  the  tellurium  antisite,  although 
to  confirm  this  possibility  a  more  quantitative  prediction  of  the  ionization  levels  of  the 
antisite  is  needed.  Deep  levels  associated  with  the  tellurium  antisite  may  also  be  respon¬ 
sible  for  the  1//  noise,  which  is  found  to  be  roughly  proportional  to  the  mercury  vacancy 
concentration®®  or  they  may  enhance  interband  tunneling  and  thereby  contribute  to  dark 
currents. 

The  role  of  native  point  defects  in  self-diffusion 

Self-diffusion  of  mercury  is  important  for  the  annealing  of  as-grown  mercury-vacancy 
doped  p— type  material  to  n— type,  to  form  p—n  junctions,*®  and  for  understanding  junc¬ 
tion  stability  in  HgCdTe  devices.  Although  both  the  mercury  vacancies  and  interstitials 
are  mobile  and  contribute  to  mercury  diffusion,®®  it  is  the  interstitial  diffusion  that  is 
found  to  dominate  in  the  modeling  of  low-temperature  anneals.*®-*®  Thus,  although  we 
find  equilibrium  mercury  interstitial  concentrations  that  are  negligible  in  terms  of  their 
contribution  to  the  net  carrier  densities,®  they  will  be  important  to  mercury  transport  in 
the  material,  and  therefore  their  properties  are  of  interest.  For  both  mercury  vacancy  and 
interstitial  diffusion,  more  than  just  defect  concentrations  enter  into  the  determination  of 
the  diffusion  coefficients;  for  the  present  discussion,  we  will  address  only  how  our  defect 
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concentrations  relate  to  the  mecisured  diffusion  coefficients. 

We  predict  that  mercury  interstitials  are  donors  and  in  equilibrium  are  present  in 
concentrations  that  cure  less  than  10®  cm~^  at  250*  C.®  We  have  compared  our  results  with 
those  of  annealing  simulations  by  the  Stanford  group^^  and  find  that  the  concentration 
of  interstitials  we  have  predicted  at  ~  200®  C  are  several  orders  of  magnitude  too  small 
to  account  for  their  modeling  of  the  formation  of  p-n  junctions.  Although  a  number 
of  approximations  in  our  calculation  of  the  formation  energies  will  affect  our  interstitial 
formation  free  energy  (for  example,  the  use  of  supercells  and  approximations  of  the  ioniza¬ 
tion  energies)  we  do  not  expect  these  to  account  for  this  large  a  discrepancy.  As  discussed 
in  the  Methodology  section  above,  we  have  completed  the  calculation  of  the  electronic 
contribution  to  the  formation  free  energies  using  pure  HgTe.  For  the  mercury  vacancy 
and  the  tellurium  antisite  that  are  the  major  defects  in  HgCdTe  and  that  occupy  the 
cation  sublattice  and  therefore  are  surrounded  by  four  tellurium  atoms,  this  is  probably 
not  such  a  bad  assumption.  Corrections  to  the  electronic  energy  due  to  the  presence  of 
cadmium  may  be  larger  for  the  interstitials  that  see  four  cations  as  their  first  nearest 
neighbors  in  one  tetrahedral  site  and  six  cations  as  near  second  neighbors  in  the  other 
tetrahedral  site.  We  are  currently  calculating  the  correction  to  our  electronic  energies, 
taking  exphcit  account  of  cadmium  in  the  lattice  to  see  if  it  will  eliminate  the  discrepancy 
with  the  Stanford  model  predictions.^* 

There  has  been  a  recent  investigation*^  in  x  =  0.22  and  0.24  Hgi_*CdrTe  on  samples 
in  which  high  concentrations  of  nonequilibrium  mercury  interstitials  have  been  introduced. 
The  observed  deep  levels  near  45  and  60  meV  above  the  valence-band  edge  were  argued 
to  be  associated  with  the  mercury  interstitials.  As  discussed  above,  at  this  point  our 
calculations  are  not  able  to  determine  the  precise  positions  of  the  defect  levels  in  the 
gap,  but  rather  have  predicted  only  that  the  mercury  interstitials  will  be  donor-like.  If  in 
fact  there  are  donor  levels  associated  with  the  mercury  interstitials  that  are  near  midgap, 
we  would  expect  them  to  exhibit  a  series  of  levels  corresponding  to  the  different  alloy 
environments  about  an  interstitial.  In  addition,  there  are  two  classes  of  interstitials  that 
we  expect  to  have  relatively  high  equilibrium  densities,  both  of  which  occupy  tetrahedral 
sites;  the  first  is  surrounded  by  four  cation  nearest  neighbors,  and  the  second  is  surrounded 
by  four  anion  nearest  neighbors.  Although  it  is  tempting  to  associate  the  two  defect  levels 
experimentally  observed  with  these  two  claisses  of  interstitial  sites,  there  is  no  apparent 
reason  why  two  levels  are  seen  in  the  x  =  0.22  material  and  only  one  level  in  for  x  =  0.24. 

Finally,  we  return  to  the  discussion  of  the  mechanism  for  tellurium  antisite  diffusion 
and  readdress  the  lockout  mechanism  that  proceeds  via  the  reaction 

Hg/  -f-  Teng  — ►  Te/.  (7) 

Our  calculations  predict  that  this  is  am  exothermic  reaction,  with  am  energy  of  ~  0.8  eV 
(the  entropy  gained  in  having  two  point  defects  rather  than  just  one  must  also  be  consid¬ 
ered  in  establishing  the  equihbrium  concentrations  of  these  defects).  In  nonequilibrium 
situations  in  which  mercury  interstitials  are  introduced  into  the  material — ^for  example, 
during  ion-beam  milling  or  oxide  bahing*® —  this  reaction  will  be  pushed  to  the  right,  and 
excess  tellurium  interstitial  will  be  produced.  Thus,  in  such  situations,  one  might  expect 
tellurium  antisite  diffusion  via  the  kickout  mechanism  to  be  a  stronger  competitor  to  the 
mercury  vacancy  mechanism,  although  the  barriers  to  the  formation  of  Hg/— Teng  pairs 
still  exist,  as  discussed  above. 

In  presenting  the  tellurium  antisite  as  a  candidate  for  the  residual  donor  we  airgued 
why  its  density  might  be  fixed  in  the  n-type  material.  On  the  other  hand,  when  discussing 
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the  possibility  that  it  is  related  to  a  Shockley-Read  recombination  center  in  vacancy-doped 
p-type  material,  we  argued  why  its  density  might  vary,  depending  on  the  cooldown  rate 
and  so  forth.  These  two  arguments  are  somewhat  inconsistent.  Although  the  possibility 
still  exists  that  the  tellurium  antisite  is  both  the  residual  donor  (via  a  first  ionization  level 
that  resonates  in  the  gap)  and  a  Shockley-Read  center  (via  a  midgap  second  ionization 
level),  to  be  convincing  a  firmer  correlation  between  the  two  would  have  to  be  established. 

The  mercury  vacancy  tellurium  antisite  pair 

Several  additional  consequences  of  the  presence  of  the  mercury  vacancy  tellurium  an¬ 
tisite  pairs  should  be  discussed.  As  one  can  see  from  Fig.  3,  we  are  predicting  a  very  large 
concentration  of  the  defect  pairs,  which  may  even  exceed  the  mercury  vacancy  concen¬ 
trations  for  LPE  material  grown  from  the  tellurium-rich  melts.  If  the  pair  is  electrically 
inactive,  as  our  preliminary  calculations  predict,  it  will  not  impact  the  carrier  concentra¬ 
tions  or  mobility.  Such  a  large  density  of  pjurs  does  imply  that  the  nonstoichiometry  of 
the  material  is  larger  than  that  due  to  the  vacancy  concentrations,  particularly  at  lower 
temperatures  where  we  predict  that  the  equilibrium  concentrations  of  neutral  pairs  will 
approach  that  of  the  mercury  vacancy. 

Note  that  the  large  binding  energy  of  the  mercury  vacancy  tellurium  antisite  pair 
suggests  that  other  bound  pairs  may  be  present  in  the  material.  For  example,  a  bound 
Frenkel  defect  (involving  the  mercury  vacancy-mercury  interstitial  pair)  that  involves 
an  acceptor  and  a  donor  defect  with  opposite  lattice  strains  may  be  important,  and  will 
impact  diffusion  of  mercury  in  the  lattice.  The  most  likely  consequences  will  be  to  increase 
the  annihilation  capture  cross  section  of  mercury  interstitials  into  mercury  vacancies  and 
to  present  a  barrier  to  the  formation  of  free  Frenkel  defects  through  a  geminate  process. 
These  phenomena  will  be  important  to  understanding  diffusion  in  HgCdTe.*^  Evidence  for 
the  mercury  vacancy  substitutional  indium  pair  has  been  seen  using  a  nuclear  hjq)erfine 
technique.*®  This  defect  is  similar  to  the  mercury  vacancy  tellurium  antisite  pair,  and  thus 
we  expect  it  may  be  well  bound.  Both  of  these  pairs  merit  further  investigation. 

DEFECTS  IN  CdTe 

Like  HgCdTe,  CdTe  has  a  wide  stability  region  and  can  be  doped  both  p-  and  n- 
type.  Its  use  as  both  a  substrate  (along  with  Cdi_,Zn,Te**)  and  a  passivant  make  it  an 
important  material  in  the  manufacturing  of  HgCdTe  LWIR  FPAs. 

We  have  predicted  the  density  of  neutred  native  point  defects  in  CdTe  as  a  function 
of  temperature  and  pressure;  results  for  700*C  are  shown  in  Fig.  4.  We  predict  the  same 
electrical  type  (donor  versus  acceptor)  for  the  native  point  defects  as  in  HgCdTe;  for 
example  we  find  that  the  cadmium  vacancy  is  an  acceptor  and  the  cadmium  interstitial 
and  tellurium  antisite  are  donors.  Based  on  the  neutral  native  point  defects  we  expect 
that  CdTe  solid  can  exist  with  both  excess  cadmium  and  excess  tellurium.  We  find  that 
the  cadmium  interstitial  is  the  dominant  defect  for  high  cadmium  pressures,  and  the 
cadmium  vacancy  is  most  important  at  mid  and  low  cadmium  pressures. 

Because  the  native  point  defects  may  have  energy  levels  in  the  gap  associated  with 
them,  it  is  important  to  include  these  in  the  analysis.  We  have  begun  to  calculate  the 
ionization  levels  associated  with  the  native  point  defects,  but  have  not  yet  included  the 
Jtihn-Teller  distortions,*®  which  can  significantly  alter  the  localized  energy  levels.  For  this 
discussion  we  shall  assume  that  the  cadmium  vacancy  has  a  first  ionization  level  near 
the  valence-band  edge  and  a  second  level  near  midgap,*®  and  assume  that  the  cadmium 
interstitial  is  a  double  shallow  donor.*  In  this  ceise  we  predict  a  p-  to  n-type  conversion 
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at  700®  C  at  relatively  high  cadmium  pressures,  in  agreement  with  experiment.^ 

From  Fig.  4  we  see  that  the  tellurium  antisite  becomes  more  important  as  the  cad¬ 
mium  pressure  is  lowered.  Depending  on  the  energies  associated  with  the  tellurium  cintisite 
donor  levels,  the  material  may  be  highly  compensated  at  the  lowest  pressures,  or  another 
p-  to  n-type  conversion  may  even  occur.  Experiments  designed  to  equilibrate  on  the 
cadmium-poor  side  of  the  stability  region  can  be  used  to  test  for  the  presence  of  tellurium 
antisites  through  both  their  impact  on  electrical  activities  and  the  presence  of  localized 
levels. 


CONCLUSIONS  AND  SUGGESTED  EXPERIMENTS 

The  theory  developed  here  clearly  indicates  that  in  Hgo.8Cdo.2Te  equilibrated  at  typ¬ 
ical  a.Tinf»a.1ing  temperatures  and  pressures  (the  right  side  of  Fig.  3(c))  the  tellurium 
antisite  and  mercury  vacancy  tellurium  antisite  pair  densities  are  well  below  levels  that 
can  impact  device  performance.  The  primary  outstanding  question  is  whether  or  not 
tellurium  antisite  diffusion  rates  are  high  enough  that  normally  processed  samples  fully 
equilibrate.  Experiments  are  needed  to  modify  these  defect  populations  in  a  controlled 
manner  so  their  impact  on  carrier  concentrations,  lifetimes,  and  other  device  properties 
can  be  determined.  The  basic  idea  is  to  modify  the  tellurium  antisite  concentration  by 
choosing  annealing  temperatures,  mercury  pressures,  and  times  that  are  long  enough  to 
permit  a  measurable  portion  of  samples  to  equilibrate. 

First  we  consider  annealing  experiments  to  test  for  the  possibility  that  the  antisite  is 
the  residual  donor.  If  tellurium  antisite  diffusion  is  so  slow  that  it  is  not  equilibrated  during 
low— temperature  mercury— saturated  anneals,  one  must  first  anneal  at  higher  temperatures 
and  lower  mercury  pressures  to  introduce  tellurium  antisite  mercury  vacancy  pairs  (this 
step  could  be  eliminated  for  LPE  material  grown  from  the  tellurium  melt),  followed  by  a 
lower-temperature,  low-mercury-pressure  anneal  to  reduce  the  antisite  density,  while  still 
maintaining  a  relatively  high  density  of  pairs,  concluded  by  a  low-temperature,  mercury- 
saturated  anti^Al  to  reduce  the  mercury  vacancy  concentration  even  further.  A  series  of 
experiments  in  which  the  conditions  of  the  first  two  anneals  were  varied  could  be  designed 
to  test  for  the  effect  they  have  on  the  residual  donor  concentration.  A  similar  experiment 
could  be  performed  to  correlate  the  Shockley-Read  center  in  mercury-vacancy-doped 
material  with  noneqmlibrium  tellurium  antisite  densities. 

Although  we  have  proposed  a  series  of  experiments  to  establish  whether  the  tellurium 
antisite  is  the  residual  donor  or  a  Shockley-Read  recombination  center  based  on  nonequi¬ 
librium  densities,  one  may  be  able  to  design  a  series  of  experiments  in  which  equilibrium 
populations  of  tellurium  antisites  are  obtained  by  choosing  high  enough  temperatures, 
thin  enough  samples,  and  long  enough  annealing  times.  One  could  choose  annealing  con¬ 
ditions  to  manipulate  the  mercury  vacancy  and  tellurium  antisite  populations  indepen¬ 
dently.  For  example  a  300®C  anneal  to  reduce  the  hole  concentration  to  ~  10*®  will  result 
in  an  order-of-magnitude  more  tellurium  antisites  than  an  anneal  at  400®  C  to  achieve  the 
same  hole  concentration;  Fig.  2  can  be  used  in  guiding  such  a  study.  The  correlation  of 
the  tellurium  antisite  densities  with  the  Shoddey-Read  center  should  be  possible  by  doing 
minority  carrier  lifetime  and/or  DLTS  measurements  on  materiads  with  the  same  mercury 
vacancy  hole  concentrations  achieved  by  anneals  at  different  temperatures.  Although  a 
study  of  this  type  was  recently  presented,*  cind  for  a  given  hole  concentration  a  corre¬ 
lation  of  lifetimes  with  the  annealing  temperature  was  observed,  the  effective  annealing 
conditions  there  were  set  by  the  stoichiometry  of  CdTe  cap  layers  and  thus  were  more 
complicated  thain  a  simple  anneal  with  a  controlled  mercury  overpressure. 

Experiments  to  test  the  presence  of  the  mercury  vacancy  tellurium  antisite  pair  are 
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similar  to  those  proposed  above  for  the  isolated  tellurium  cintisite.  Our  preliminary  cal¬ 
culations  indicate  that  the  pair  will  be  electrically  inactive  and  have  no  states  in  the  band 
gap.  As  such,  we  do  not  expect  their  presence  to  have  a  direct  impact  on  the  electrical 
properties,  but  their  presence  will  be  manifested  in  their  impact  on  tellurium  antisite 
diffusion.  Note  that  these  experiments  proposed  to  test  for  the  presence  of  the  tellurium 
antisite  cannot  discriminate  between  the  isolated  tellurium  antisite  and  the  mercury  va¬ 
cancy  tellurium  antisite  pair. 
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FIGURE  CAPTIONS 


Fig.  1.  Hole  concentrations  at  77  K  as  a  function  of  mercury  partial  pressure  for 
material  annealed  at  high  temperatures,  as  indicated.  Experimental  data  were  taken 
from  Vydyanath.®  Results  of  our  theory  are  shown  as  solid  lines.  To  demonstrate 
the  sensitivity  of  our  predictions  to  small  changes  in  the  calculated  formation  free 
energies,  we  show  the  dashed  lines,  which  are  the  theoretical  results,  but  with  the 
mercury  vacancy  formation  energy  increased  by  10%  tind  including  a  rigid  upward 
shift  of  the  hole  concentrations  by  a  factor  of  5.5.  Note  that  the  results  shown  here 
differ  from  those  in  Ref.  9  because  the  inclusion  of  the  additional  degeneracy  factor 
for  the  singly  ionized  tellurium  antisite. 

Fig.  2.  Total  mercury  vacancy  (solid  line)  and  tellurium  antisite  (dashed  line)  densi¬ 
ties  as  a  function  of  mercury  partial  pressures  (atm).  Full  equilibration  of  all  defects 
is  assumed  at  all  temperatures  (see  text  for  further  discussion). 

Fig.  3.  The  concentration  of  native  points  as  a  function  of  mercury  pressure  within 
the  phase  stability  region  at  (a)  500®C,  corresponding  to  the  LPE  growth  temperature, 
(b)  185®C,  the  growth  temperature  for  MBE,  and  (c)  220®C,  a  typical  temperature 
for  mercury-saturated  anneals.  A  range  of  concentrations  for  the  mercury  vacancy 
tellurium  antisite  pair  is  shown,  based  on  our  preliminary  results,  as  discussed  in  the 
text. 

Fig.  4.  Neutral  native  point  defect  densities  of  CdTe  within  the  stability  region  at 
700®C. 
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Abstract 


Band-edge  shifts  induced  by  the  electron-phonon  interaction  are  calcu¬ 
lated  for  HgCdTe  alloys  and  various  semiconductor  compounds  starting  from 
accurate  zero-temperature  band  structures.  The  calculated  temperature  vari¬ 
ation  of  gaps  agrees  with  experiments  to  better  than  10%  in  all  materials 
except  InAs  and  InSb  where  the  deviation  is  about  50%.  While  the  simple 
picture  that  the  intra  (inter)-band  transitions  reduce  (increase)  the  gap  still 
holds,  we  show  that  both  the  conduction  band  edge  Ec  and  valence  band  edge 
Ev  move  down  in  energy.  These  shifts  in  Ey  affect  the  valence  band  offsets  in 
heterojunctions  at  finite  temperature.  The  temperature  variations  of  valence 
band  offset  and  the  electron  effective  mass  are  also  reported. 

Key  words:  Temperature-dependent  band  structures,  electron-phonon  in¬ 
teractions,  band  offset,  HgCdTe  and  alloys,  III-V  semiconductors 
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INTRODUCTION 


The  temperature  (T)  dependence  of  energy  gaps  of  semiconductors  is  of  great  physical 
and  technological  interest.  The  quantities  such  as  band  offset  and  effective  mass  depend 
sensitively  on  the  temperature  variation  of  band  edges.  Numerous  theoretical  [1-9]  and  ex¬ 
perimental  [10-19]  studies  have  been  undertaken  to  obtain  both  qualitative  and  quantitative 
variations  of  various  gaps  in  semiconductors.  The  gap  decreases  with  increasing  tempera¬ 
ture  in  medium-gap  and  wide-gap  semiconductors,  and  it  increases  in  small-gap  materials 
such  as  HgCdTe,  PbS,  PbSe,  and  PbTe.  The  thermal  expansion  of  the  lattice  and  electron- 
phonon  interactions  are  usually  considered  causes  for  the  temperature  variation  of  the  band 
structures.  Thermal  expansion  always  reduces  gaps. 

In  a  perturbation- theory  treatment  of  electron-phonon  interactions,  the  intraband  tran¬ 
sitions  reduce  the  gap  whereas  interband  transitions  increase  it,  and  the  net  shift  in  the 
gap  ran  be  positive  or  negative.  Here  we  calculate  the  gap  variation  in  Hgi_iCdj;Te  alloys, 
GaAs,  InAs,  InSb,  InP,  and  CdTe  compounds,  starting  from  accurate  band  structures,  wave 
functions,  proper  phonon  dispersion  relations,  and  taking  account  of  matrix  elements  of  the 
electron- phonon  interactions.  The  contributions  from  each  phonon  branch  to  each  electron 
band  have  been  obtained  to  assist  physical  understanding  of  the  underlying  causes  of  the 
variations.  We  show  that  both  conduction  and  valence  band  edges  move  down  in  energy. 
When  the  valence  band  edge  moves  more  than  the  conduction  band  edge,  the  gap  increases 
with  T,  as  in  the  case  of  some  Hgi_*Cd3,Te  alloys  with  x  <  0.5.  The  reverse  occurs  for  all 
III-V  compounds  studied  and  Hgi-xCd^Te  with  x  >  0.5.  This  observation  has  an  important 
effect  on  our  understanding  of  the  variation  of  band  offsets  in  semiconductor  heterojunc¬ 
tions.  In  addition  to  the  gap,  other  features  of  the  band  structure  change  with  temperature 
and  will  affect  the  spectral  variations  of  the  absorption  coefficient  and  transport  properties. 
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METHOD 


Our  calculation  of  the  temperature  dependence  of  the  band  gap  starts  with  accurate  band 
structures.  Empirical  pseudopotential  form  factors  are  used  to  construct  a  hybrid  pseudopo¬ 
tential  tight-binding  (HPTB)  Hamiltonian.  The  pseudopotential  part  of  the  Hamiltonian  is 
universal-  it  applies  to  all  group  IV,  III-V  and  II- VI  compounds.  The  smaller  tight-binding 
part  is  expressed  in  a  minimum  set  of  sp^  Slater  orbitals  per  atom.  This  Hamiltonian  is 
then  transformed  into  an  orthonormal  basis.  A  site-diagonal  spin-orbit  Hamiltonian  is  then 
added.  Parameters  in  the  tight-binding  perturbation  are  chosen  to  fine-tune  the  band  struc¬ 
tures  to  agree  well  with  experiments  [20,21].  Various  results  obtained  using  these  band 
structures  are  found  to  be  quite  reliable  [21-23].  The  present  study  subjects  the  accuracy  of 
the  wave  functions  as  well  as  the  energies  to  a  sensitive  test. 

The  dilation  contribution  to  the  band  gap  reduction  is  given  [5,9]  by  ZarBOEgld?, 
where  the  thermal  expansion  coefficient  of  the  lattice  aj,  the  bulk  modulus  B,  and  the 
change  in  the  gap  with  pressure  are  obtained  from  the  hterature  [19].  The  electron- phonon 
interactions  with  all  phonon  branches  that  cause  the  band  structure  changes  are  treated 
in  perturbation  theory.  The  total  Hamiltonian  is  assumed  to  be  a  sum  of  potentials  from 
single  atoms.  The  atomic  potential  in  the  solid  is  traditionally  expanded  in  a  Taylor  series, 
with  only  the  leading  term  retained,  and  the  energy  shifts  it  causes  are  evaluated  in  second- 
order  perturbation  theory.  However,  it  has  been  demonstrated  by  a  number  of  researchers 
[3,4,6]  that  retention  of  first-order  perturbation  terms  with  a  second  term  in  the  Taylor 
series  expansion  is  necessary  to  preserve  symmetry.  We  retain  both  terms.  The  change  in 
the  energy  at  a  given  wave  vector  k  is 


n'k' 


<n^k1Vi|nk> 
£^nk  -  -En'k' 


(1) 


where  Vi  and  V2  are  the  first  two  terms  in  the  Taylor  expansion  of  the  total  electron-phonon 
potential  in  powers  of  atomic  displacements  In  the  TB  formalism,  Eq.  (1)  can  be  written 
in  terms  of  the  matrix  elements 
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<  \'j'c,'\Vi\\ja  >=  VV,„.(dii!')  ■  (fly  -  fli), 


(2) 


and 

<  I'jWmija  >  =  l[f,.y  •  (V)  Vo^.Cd;/:')  ■  fiy 

+fir(V)V„.(d;i':')-f,;],  (3) 

where  d/ji*  is  the  position  vector  connecting  atomic  sites  1,  species  (anion  or  cation)  j  and  site 
1',  species  j',  and  Yaa'  {^(v)  is  a  HPTB  matrix  element  between  the  orbitals  a  and  a'  located 
on  those  atoms.  From  the  quantum  theory  of  harmonic  crystals,  the  atomic  displacements 
^  can  be  expressed  in  terms  of  normal  modes-  that  is,  phonons.  We  have 

(4) 

where  q  and  w  are  phonon  wave  vector  and  frequency,  A  denotes  phonon  branch,  a(al)  is 
a  phonon  annihilation  (creation)  operator,  M  is  the  atomic  mass,  and  e  is  an  eigenvector 
in  a  diamond  or  zincblende  structure  of  the  six-dimensional  dynamical  matrix  eigenvalue 
problem 

Mu?e  =  D(q)e  (5) 

Evaluation  of  the  matrix  elements  given  by  Eqs.  (2)  and  (3)  requires  knowledge  of  spatial 
variations  of  the  interatomic  TB  matrix  elements  Yaa'-  Iii  Harrison’s  universal  TB  approach 
[24],  these  matrix  elements  scale  as  d"^.  In  our  generalization,  we  assume  that  Yaa'  varies 
as  d"*"  and  the  repulsive  first-neighbor  pair  energy,  following  Harrison’s  overlap  argument, 
as  7;/d^’".  The  two  unknowns  m  and  7/  are  determined  by  requiring  that  the  calculated 
equilibrium  bond  length  and  bulk  modulus  agree  well  with  experiments.  This  approach, 
with  electrons  and  phonons  treated  from  the  same  underlying  Hamiltonian,  has  previously 
been  used  successfully  to  explain  hot  electron  transistor  characteristics  [23]  and  is  also  in 
fairly  good  agreement  with  first-principles  calculations  [25].  The  dynamical  matrix  D  is 
calculated  from  the  valence  force  field  model  [26]. 
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The  calculational  procedure  is  as  follows.  For  a  chosen  material,  m  and  ??  are  evaluated. 
Then  the  first  and  second  derivatives  of  all  interatomic  matrix  elements  are  obtained.  The 
dynamical  matrix  is  diagonalized  to  obtain  w  and  e  as  a  function  of  q  and  A.  The  phonon 
structures  and  electronic  band  structures  are  used  [Eqs.  (1)  through  (4)]  to  obtain  the  change 
in  the  band  energy  at  a  given  k.  The  polar  coupling  terms  are  included  in  the  longitudinal 
optical  phonon  contributions.  When  we  are  interested  in  studying  the  change  in  the  direct 
gap,  k  is  taken  to  be  zero.  However,  when  the  temperature  variation  of  the  effective  masses 
or  indirect  gap  are  studied,  non-zero  k  values  must  be  used  and  the  Brillouin  zone  summation 
in  Eq.  (1)  should  be  carried  over  the  entire  zone  with  reduced,  or  no,  symmetry. 


RESULTS 

The  calculated  band-gap  change  as  a  function  of  T  in  Hgo.r8Cdo.22Te  is  shown  in  Fig. 
1,  With  increasing  T,  the  direct  gap  increases  in  Hgo.78Cdo.22Te.  Notice  that  the  calculated 
values  are  typically  within  10  to  15  meV  of  experimental  values  [16-18].  The  cross  (x)  at 
T=0  represents  the  calculated  zero-point  correction  to  the  gap  (13.6  meV  for  H5^o.78Cdo.22Te). 
The  zero-temperature  band  gap  calculated  without  electron-phonon  interactions  should  have 
this  correction  subtracted  for  comparison  to  experimental  values. 

The  change  in  the  gap  is  traditionally  explained  in  terms  of  inter-  and  intraband  interac¬ 
tions.  The  intravalence  (conduction)  band  interactions  push  the  valence  (conduction)  band 
edge  up  (down),  thus  reducing  the  gap.  Similarly,  the  valence-conduction  band  interactions 
increase  the  gap.  Hence,  one  might  expect  the  gap  to  decrease  in  wide-gap  semiconductors 
and  possibly  increase  in  small-gap  semiconductors.  In  addition,  arguments  based  only  on 
total  density  of  states  and  ignoring  variations  in  matrix  elements  will  predict  the  valence 
band  edge  Ey  move  up  in  energy,  because  the  hole  effective  mass  is  one  to  two  orders  of 
magnitude  larger  than  the  electron  mass.  As  seen  from  Fig.  2,  our  detailed  calculations  of 
band  edge  movements  in  Hgo.78Cdo.22Te  do  not  support  this  traditional  view.  We  find  that 
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at  high  temperatures  (>150  K).  The  calculated  dEj/dT  values  (circle)  are  compared  with 
experiments  (cross)  values  in  Fig.  3.  We  see  that  the  calculations  produced  correct  trends  in 
all  these  materials,  but  compare  less  favorably  with  experiments  in  InAs  and  InSb.  However, 
it  is  important  to  note  that  the  sign  of  the  change  is  not  exclusively  determined  by  the 
magnitude  of  the  zero-temperature  gap.  For  example,  although  Hgo.70Cdo.30Te  and  InSb 
have  the  same  zero-temperature  gap  of  0.235  eV,  the  InSb  gap  decreases  with  T,  whereas 
the  Hgo.7oCdo.3oTe  gap  increases  with  T.  The  combination  of  gap  size,  conduction  band 
width,  and  intervalley  separations  gives  rise  to  these  interesting  variations  in  the  gap  with 
T. 

The  observation  that  both  Ey  and  Ec  move  down  in  energy  has  an  important  effect  on 
band  offsets  in  heterojunction-based  devices.  For  example,  the  zero-temperature  valence 
band  offset  between  Hgo.78Cdo.22Te  and  CdTe  is  believed  to  be  around  350  meV.  However, 
we  find  that  at  300  K,  Ey  in  Hgo.78Cdo.22Te  and  in  CdTe  moves  down  by  215  meV  and 
30  meV,  respectively.  If  the  dipole  contribution  remains  the  same,  the  valence  band  off¬ 
set  decreases  to  165  meV  at  300  K.  The  contention  that  the  dipole  contribution  is  nearly 
temperature  independent  stems  from  the  observation  that  any  shift  in  the  average  effective 
crystal  potentialshould  effectively  be  screened  out,  since  these  semiconductors  are  good  di¬ 
electrics  (c  >  10).  The  temperature  variation  of  the  bands  should  be  taken  with  respect  to  a 
fixed  average  potential.  For  our  Hamiltonian,  the  valence  band  edge  movements  in  each  side 
of  the  junction  are  calculated  with  reference  to  a  fixed  average  state.  Thus,  the  calculated 
temperature  dependence  of  the  difference  in  the  VB  edge  of  the  constituent  hetero junction 
materials  effectively  governs  the  temperature  dependence  of  the  band  offset.  In  addition  to 
the  electron-phonon  interactions  discussed  above,  lattice  dilation  changes  the  band  edges 
differently  [9].  This  effect  is  not  included  here.  In  any  case,  this  band  offset  change  has 
important  imphcations  for  the  design  of  abrupt  heterojunction  IR  absorption  and  confined 
well  laser  devices. 

In  principle,  the  band  structure  at  any  wave  vector  k  will  change  with  temperature.  With 
the  change  in  the  fundamental  gap,  the  band  curvature  (or  effective  mass)  also  changes  thus 
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affecting  the  optical  and  transport  properties  of  the  material.  The  self-energy  calculated 
in  this  method  will  include  the  effect  of  scattering  due  to  phonons  and  the  change  in 
the  temperature-dependent  band  structure  self-consistently.  In  the  case  a  of  narrow-gap 
material  such  as  Hgo,78Cdo.22Te,  the  effective  mass  alone  does  not  explain  the  low-energy 
portion  of  the  conduction  band  structure.  The  lowest  CB  energy  at  any  k  is  best  described 
by  a  hyperbola  [27],  (7k^  +  -  c.  The  calculated  band  gap,  effective  mass,  7,  and  c  as 

functions  of  T  are  given  in  Table  2.  The  effective  mass  and  c  are  directly  proportional  to 
the  gap  and  hence  monotonically  increase  with  T.  This  is  expected  from  the  k  •  P  theory 
argument,  but  the  magnitudes  predicted  by  two  theories  differ.  7  decreases  slightly  at  lower 
temperatures  and  then  starts  to  increase  with  T.  In  a  previous  publication  [27],  we  had 
simulated  these  temperature  variations  of  7  and  c  by  adjusting  the  Hg  concentration  in 
HgCdTe  alloys  to  produce  proper  gap  at  each  temperature.  Those  values  are  in  remarkable 
agreement  with  the  values  reported  in  Table  2.  We  conclude  that  7  and  c  (given  in  Table 
2)  can  be  interpolated  to  considerable  accuracy  for  any  positive  gap  in  the  HgCdTe  alloys 

CONCLUSIONS 

Although  the  calculations  produced  correct  trends  in  all  materials,  the  calculated  changes 
in  the  band  gap  of  InAs  and  InSb  were  about  a  factor  of  2  smaller  than  in  the  experiments. 
We  find  that  our  calculated  TA  phonon  frequencies  away  from  zone  center  in  these  com¬ 
pounds  were  considerably  larger  than  those  found  in  experiments.  As  noted  from  Table 
1,  a  substantial  contribution  comes  from  acoustic  phonons.  Consequently,  our  theoreti¬ 
cal  values  of  Ej(T)  are  smaller  than  in  experiments.  Better  predictability  should  result 
from  improvement  in  the  dynamical  matrix  calculated  from  the  underlying  Hamiltonian.  In 
addition,  at  higher  temperatures  higher-order  perturbation  terms  must  be  included  along 
with  finite-temperature  ‘renormalized’  bands  rather  than  the  zero-temperature  bands.  Such 
renormalization  affects  the  monotonic  change  in  the  gap  and  introduces  nonlinear  terms. 

In  summary,  we  have  calculated  the  temperature  variations  of  band  gaps  in  various  semi- 
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conductors.  A  fairly  accurate  HPTB  Hamiltonian  is  used  in  the  calculation  of  electron  and 
phonon  structures.  The  calculations  explain  the  increase  in  the  band  gap  of  Hgo.js  Cdo.jjTe, 
and  the  decrease  in  the  band  gap  of  all  IITV  compounds  studied.  We  show  that  acoustic 
phonons  make  the  major  contribution.  Contrary  to  traditional  thinking  based  on  total  den¬ 
sity  of  states  arguments,  we  find  that  both  the  valence  and  the  conduction  band  edges  move 
down  in  energy.  One  important  consequence  of  this  observation  will  be  m  the  band  offsets 
in  semiconductor  hetero junction  devices.  Finally,  there  is  a  small  and  usually  negligible 
zero-point  motion  contribution  to  low-temperature  band  gaps  arising  from  electron-phonon 
interactions. 

We  thank  Dr.  M.  Cardona  of  the  Max  Planck  Institute,  Stuttgart,  for  pointing  us  to 
several  references.  Funding  from  ONR  (contract  N00014-93-C-0091)  and  ARPA  (contract 
MDA972-92-C-0053)  is  gratefully  acknowledged. 
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TABLES 


TABLE  1.  Calculated  change  in  the  valence  (v)  and  conduction  (c)  band  edge  energies  (in 
meV)  of  Hgo.78Cdo.22Te  alloy.  Contributions  from  interaction  with  various  phonon  modes 
are  shown  in  rows  3  to  10. 


band 

1 

2 

3 

4 

5 

6 

7 

8 

V 

5.90 

13.99 

56.85 

88.51 

-80.26 

-92.49 

-97.45 

-102.07 

Total 

c 

1.04 

2.84 

17.27 

34.07 

-23.81 

-25.43 

-43.96 

-42.48 

V 

4.03 

6.21 

30.06 

41.52 

-49.86 

-62.17 

-61.04 

-49.62 

TA 

c 

0.26 

0.65 

12.41 

27.47 

-3.58 

-14.81 

-24.50 

-15.40 

V 

0.91 

3.18 

3.14 

10.53 

-10.63 

-5.92 

-11.41 

-20.59 

LA 

c 

0.56 

1.29 

1.27 

1.72 

-14.94 

-3.00 

-2.66 

-7.17 

V 

0.33 

2.11 

6.12 

11.22 

-9.52 

-4.98 

-11.80 

-17.40 

LO 

c 

0.13 

0.63 

0.75 

0.44 

-3.78 

-2.45 

-1.18 

-7.31 

V 

0.63 

2.48 

17.53 

25.24 

-10.26 

-19.42 

-13.20 

-14.45 

TO 

c 

0.10 

0.27 

2.84 

4.44 

-1.52 

-5.16 

-15.62 

-12.59 
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TABLE  2.  Calculated  Eg  [meV],  7  [eV^],  c[eV]  and  the  effective  mass  ratio  of  Hgo.78Cdo.22Te 
alloy  as  functions  of  temperature.  m*(0)  is  0.008.  The  zero  point  correction  is  13.6  meV. 


T 

7 

c 

m*(T)/m*(0) 

1.00 

113.60 

47.7656 

0.0588 

1.0000 

10.00 

112.67 

47.7553 

0.0588 

1.0005 

20.00 

112.56 

47.7169 

0.0592 

1.0072 

30.00 

114.44 

47.6421 

0.0598 

1.0199 

40.00 

117.15 

47.5582 

0.0607 

1.0361 

50.00 

120.42 

47.4461 

0.0615 

1.0532 

60.00 

123.96 

47.3310 

0.0624 

1.0714 

70.00 

127.65 

47.2091 

0.0634 

1.0904 

80.00 

131.44 

47.0821 

0.0643 

1.1095 

90.00 

135.28 

46.9418 

0.0653 

1.1288 

100.00 

139.17 

46.7964 

0.0662 

1.1483 

150.00 

158.85 

46.1544 

0.0712 

1.2529 

200.00 

178.73 

45.5930 

0.0767 

1.3669 

250.00 

198.68 

45.2441 

0.0832 

1.4938 

300.00 

218.66 

45.1167 

0.0908 

1.6342 

350.00 

238.65 

45.3460 

0.1000 

1.7913 

400.00 

258.66 

46.0193 

0.1115 

1.9672 

450.00 

278.67 

47.3751 

0.1263 

2.1657 

500.00 

298.69 

49.8338 

0.1468 

2.3919 

550.00 

318.71 

54.0581 

0.1763 

2.6491 

600.00 

338.74 

61.7125 

0.2238 

2.9445 
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FIGURES 


FIG.  1.  Change  in  the  band  gap  of  Hgo.78Cdo,22Te  with  temperature. 

FIG.  2.  Variation  of  conduction  (dashed  line)  and  valence  (solid  line)  band  edges  of 
Hgo.78Cdo.22Te  with  temperature. 

FIG.  3.  Derivative  of  direct  gap  with  temperature  for  various  semiconductor  compounds  and 
alloys  as  a  function  of  zero-temperature  gap.  The  vertical  lines  represent  Hgo.78Cdo.22Te,  InSb, 
InAs,  Hgo.5Cdo.5Te,  InP,  GaAs  and  CdTe,  respectively. 
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